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Abstract. In the vicinity of a massive object of various scales (ranging from 
young stars to galactic nuclei) , mass flow creates a spectacular structure combining 
a thin disk and coUimated jet. Despite a wide range of scaling parameters (such as 
Reynolds number, Lundquist number, ionization fractions, Lorentz factor, etc.), 
they exhibit a remarkable similarity that must be dictated by a universal principle. 
A generalized Beltrami condition has been formulated as a succinct representation 
of such a principle. The singularity at the center of the Keplerian rotation forces 
the flow to align with the "generalized vorticity" (including the effect of localized 
density and finite dissipation) which appears as an axle penetrating the disk, 
i.e. the jet is a Beltrami flow. Based on the Beltrami flow model, an analytical 
expression of a disk-jet system has been constructed by the method of similarity 
solution. 
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1. Introduction 

The combination of a thin disk and coUimated jet is a common structure that is 
created in the vicinity of a massive object O [1] [121 E US]. Beneath a large variety of 
scales, constituents, and local processes of such systems, there must be a simple and 
universal principle that dictates the remarkably similar geometry; see Fig.[TJ Here we 
show that the coUimated structure of jet is a natural consequence of the alignment 
of the flow velocity and the vorticity, i.e. so-called Beltrami condition determines the 
structure. On a Keplcrian thin disk, the vorticity becomes a vertical vector with a 
magnitude oc r~'^/^ (r is the radius from the center of the disk), which appears as a 
spindle of the disk. Then, the alignment is the unique solution for avoiding singularity 
of Coriolis force near the center. However, we need to generalize the "vorticity" to 
deal with the strong heterogeneity of the disk-jet system, as well as to account for 
the dissipation that causes accretion. The mission of this study is to formulate an 
appropriate generalized vorticity to which the disk-jet flow aligns. 

Let us start by a short review of the Beltrami vector field which has wide 
applications as a model of various vortex structures found in nature. The Beltrami 
condition, demanding the alignment of flow and its vorticity, forces the total energy 
density (consisting of thermal energy, kinetic energy, and other energies of coupled 
fields such as gravitational or electromagnetic) to distribute homogeneously (so called 
Bernoulli condition); the Beltrami-Bcrnoulli condition, thus, fits the notion of "relaxed 
state." While a Beltrami field was discovered by many researchers as a particular 
type of equilibrium state [4l 114). or free decay solution |10|. in fluids or plasmas, its 
relation to the helicity was noticed in the study of "force-free magnetic flelds" in 
plasmas; Woltjer |27j invoked the helicity as a constraint in minimizing the magnetic 
energy; Taylor |25| considered that the "relaxed state" is the energy minimizcr under 
the lugged constancy of the total helicity; the corresponding Euler-Lagrange equation 
becomes the eigenvalue problem of the curl operator; see [29| for the mathematical 
characterization of the curl operator and its eigenfunctions. We can actually observe 
Taylor relaxed states on various experiments |26j , as well as in some astronomical 
systems; e.g. [HI [11]. The helicity constraint causes a finite vorticity in the relaxed 
state, resulting in interesting topological properties of field-lines [T^- In the context 
of Hamiltonian mechanics, the helicity is regarded as a Casimir element representing 
the defect of the governing symplectic geometry |18j ; a Beltrami field is an equilibrium 
point on a helicity leaf. 

The Beltrami fields constitute an interesting, widely-applied class of vectors (or 
axial vectors) with "twisted" field- lines; e.g. [3J [51 [53]. Not only for describing 
equilibrium states, they are applied to analyze waves [501 17] . instabilities |31|. and 
turbulence [20l [9] [28] (here we can cite only a short list of references) . 

A variety of generalizations have been proposed. Including the cross helicity 
as an additional constraint, we obtain a fiow parallel to the magnetic field and the 
resultant hydrodynamic pressure balancing with a static pressure |23| . The two-fluid 
(Hall MHD) formulation elucidates a fundamental structure in the coupling of flow 
and magnetic flcld in terms of the canonical vorticities |15j : the simultaneous ion 
and electron Beltrami conditions yields the double Beltrami fields, which have various 
applications in both laboratory and astrophysical plasmas; e.g. [El [TBI [32] . Another 
generalization is made by including the coupling of vortex and compressible motions; 
by boosting such generalized Beltrami fields, we obtain modulating nonlinear Alfven 
waves [33] . 
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Figure 1. In a disk-jet system, tiie accreting flow and jet align parallel to a 
generalized vorticity. This figure shows the streamlines of a "generalized Beltrami 
fiow" to be constructed in Sec. [5] 



This work extends the scope of Beltrami fields to show that the disk-jet system 
is a "generalized" Beltrami vortex [3T] ; the generalization is made by introducing a 
new generalized vorticity that combines the vorticity of a "reduced momentum" — the 
reduction is to account for a viscous dissipation as well as to subtract the centrifugal 
force of the Keplerian rotation. 

2. Generalized vorticity and Beltrami condition 

To demonstrate how the alignment condition arises and how it determines the 
structure of a thin disk and narrowly-coUimatcd jet, we invoke a simple model of 
neutral fluid; for a generalization to magneto-fluid, see Appendix A. Let P = pV 
denote the mechanical momentum density, where p is the mass density and V is the 
(ion) flow velocity. The momentum equation read as 

dtP + v ■ {VP) = -pV(j) - Vp - V • n, (1) 

where </> is the gravity potential, p is the scalar pressures, and 11 is the (effective) 
viscosity tensor. The variables are normalized as follows: We choose a representative 
flow velocity Vq and a mass density pQ in the disk, and normalize V and p by these 
units. The energy densities p(f) (gravitational) and p (thermal) are normalized by the 
unit kinetic energy density £o := poVq/2. The independent variables (coordinate x 
and time t) arc normalized by the system size Lq and the corresponding transit time 
To = Lo/Vo. 

We consider stationary solutions; putting 9t = in ([T]), we obtain 

V • (VP) = -p\7(j) - - V • n. (2) 
In order to derive a term that balances with the viscosity term, we decompose the 
"inertia term" [the left-hand side of as follows: we first write 

P = PlP2, (3) 
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Pi:=piV, P2:=P2V. (4) 
Using these variables, we may write 

V ■ (VP) = V ■ {pVV) = (V ■ Pi)P2 + (Pi • V)P2. (5) 

In the conventional formulation of fluid mechanics, we choose /O2 = 1 and pi = p. 



Then, using the mass conservation law dtp + V ■ P = 0, we can rewrite the left-hand 
side of (H)) as p[dtV + {V ■ V)V]. In the preset analysis, however, we choose a different 
separation of the inertia term in order to match (V • Pi)P2 with the viscosity term. 
By a reduced p2 (< 1), we will define a generalized vorticity of a reduced momentum. 
Multiplying (/52/pi) on both sides of ([2]), we obtain (assuming a barotropic relation, 
we put Vp = pV/i with an enthalpy h) 

P2 X f22 = \ypi + pIv {(p + h) + ^ [(V • Pi)P2 + V • n], (6) 

2 pi 

where 

rj2 := V X P2 (7) 

is a generalized vorticity. 

In the next section, we will show that a generalized Beltrami condition, demanding 
that $72 parallels P2 (thus, ^ holds with both sides being zero), is a unique recourse 
to avoid singular energy densities in a disk-jet geometry. 

3. Beltrami Model of Disk-Jet System 

Now we consider an axisymmetric (de = in the r-9-z coordinates) disk-jet system. 
A massive central object produces (j> — ~MG/r (we neglect the mass in the disk 
and jet). In the disk, V « Veeg with the Keplerian velocity Vs oc r~^/^. Then, 
V X V oc r^'^/^Ej,. The momentum is strongly localized in the thin disk, and the 
vorticity diverges near the axis. This particular configuration poses strong constraints 
on the force balance equation ([S]), allowing only a special class of solutions to exist; 
following conclusions are readily deducible. 

3.1. Balance of viscosity force and partial inertia 

In the disk, a radial flow (much smaller than Vo) is caused by the viscosity (a finite 
dissipation breaks the conservation of the angular momentum and enables the flow to 
cause accretion). 

Since the flow Y is primarily in the azimuthal {0) direction, the viscosity force 
can be approximated as (under assumption of the azimuthal symmetry, V • V ~ 0) 

- V • n w -V X (p?7V X y), (8) 

where 77 is the shear viscosity coefficient. In a Keplerian thin disk, where V ~ 
V(3f "^^^ee, and \/{pri) is approximately vertical, we may estimate 

- V ■ n w -pi]V X (V X y) = -pT]V„^r-^/^eg. (9) 

Hence, we may write —V • 11 = —uP with a positive coefficient i^(r), i.e., the viscosity 
force is primarily in the azimuthal (toroidal) direction, which can be balanced by the 
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term {P2/P1) (V ■ Pi)P2 that has been extracted from the inertia term; see ([5]). Using 
the steady-state mass conservation law V • P = 0, we observe 

pr'V ■ Pi = P2V ■ \/p^^ = V • Vlogp2. 

Hence, the balance of the viscosity and the partial inertia term demands 

V ■V\ogp2=v, (10) 

which determines the parameter p2. (Since dg = 0, we can integrate (jlOp for p2 along 
the streamline of V on the poloidal r-z plane.) 

The remaining part pi of the density is determined by the mass conservation law: 
By V ■ P = V ■ (pV) = P2V ■ Vpi + piV ■ Vp2 + /O1/O2V • y = and we obtain a 
relation 

V ■\7\ogpi = -\7 -V -1^. (11) 



3.2. Beltrami condition 

Near the axis, the poloidal component of the flow begins to have an appreciable vertical 
(z) component — this is the place where the jet is created; we are going to unearth 
the mechanism that collimates the flow. 

After balancing the third and fourth terms in the remaining terms do not 
have a toroidal (azimuthal) component. In fact, the right-hand-side gradient terms 
have only poloidal components, and hence, the left-hand-side P2 x $12 must not have 
a toroidal component (to put it in another way, wc have extracted the partial inertia 
term —{p2/pi){^ ■ Pi)P2 from the total inertia to separate the toroidal component). 
The vorticity Q,2 includes a singular factor V x V oc r~^/'^ez. To eliminate the 
divergence of P2 x ^2 near the axis, P2 must align to $72, i-e., the Beltrami condition 

^2 = XP2 (12) 

must be satisfied, where A is a certain scalar function. The flow V = P2/P2 is, 
therefore, collimated by the generalized vorticity r22 creating a jet. 

Remark 1 Here the essential part of the Beltrami condition is its poloidal component, 
which dictates the poloidal flow so as to eliminate the toroidal component of the 
inertia term P2 x $12- As for the toroidal flow Vgeg, which yields primarily a radial 
(centrifugal) inertia force, the "Beltrami condition" brings about an extra constraint. 
In fact, if we were to use the conventional vorticity Jl = V x V (i.e., if p2 ~ Ij 
and estimate the centrifugal force of the Keplerian flow V = Vor~^^^eg , the term 
P X ft contributes a half of the total centrifugal force (V • V)P = pV^ /r, while the 
term on the right-hand side of ^ contributes the remaining half; combining 

these two terms, we obtain the right balance with the gravity —pMG/r^. Hence, the 
conventional Beltrami condition P xfl ~ would lead to an inadequate estimate of the 
toroidal flow. However, our generalized Beltrami condition, based on the generalized 
vorticity f22 ( including P2 7^ 1 A can be made consistent with the Keplerian velocity; 
see 5'ec.O 



3.3. Bernoulli condition 

When the Beltrami condition eliminates the left-hand side of ([S]), the remaining 
potential forces must balance to achieve the Bernoulli condition [15| that reads as 

^VP^ + S/{(f) + h) = V f^V^ + 4, + h) +y2viogp2=0. (13) 
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The system of determining equations is summarized as follows: By (|10p . we 
determine p2 for a given ly. This equation involves V = P/p that is governed 
by the Beltrami equation After determining V and p2, we can solve the 

Bernoulli equation (|13p to determine /i (the gravitational potential is approximated 
by = -7\fG'/r). 



4. Parameterization by Clebsch potential 

We may rewrite the determining equations (|10p - (jl3p in a succinct form by invoking the 
Clebsch parameterization. In an axisymmetric geometry, the divergence-free vector P 
may be parameterized as 

P = Vi' X V9 + IWe, (14) 

where / = prVe- Both ijj and / do not depend on 9. Since P ■ Vi/' = 0, the level sets 
(contours) of -0 are the streamlines of P (or those oi V = P / p). In a disk region, 
rVg oc r^/^, while p is a strongly localized function with respect to z. 
Substituting ([H]) into dTU]) yields 

V =-P ■V\ogp2 = —{\ogp2,i'], (15) 
p rp 

where {a, 6} := {drb){dza) — {dra){dzb). For a given set of P, p and ly, we can solve 
(jlSp to determine p2, as well as pi ^ p/p2 that is consistent to (fTTj) . 

Substituting (|14p into the Beltrami condition (|12p . we obtain, from the toroidal 
component, 

Ap^iVV - V {p^'l) , (16) 

implying that p^^I =: I2 = hiip) and Xp^^ — /2(V') (we denote f'{ip) — df{tp)/di{j), 
and, from the poloidal component, 

/:v-vv-viogpi = -p?/^(V')/2(V'), (17) 

where Cip := rdr{r~^dr'4') + (^zV'- This elliptic partial differential equation determines 
the poloidal- momentum Clebsch potential ip. 

The Beltrami condition has decoupled the gradient forces from the momentum 
equation, which must balance separately — the Bernoulli condition (1131) which now 
reads as the determining equation of the enthalpy: 

v/i = -vf:i(|V'|' + /'^ ^ 



i'\'+I')\7\ogp2. (18) 



5. Analytic Similarity Solution 

5.1. A similarity solution modeling disk-jet structure 

In this section, we construct a similarity solution of the model (|17p . which describes 
a fundamental disk-jet structure. We define 

T := - (r > 0), (19) 
r 

and an orthogonal variable (Vt • Vcr = 0) 



a := Vr^ + z^ • (20) 
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Figure 2. The momentum field (contours of i/' that describe the streamlines 
of the poloidal component of P) of the similarity solution (with D = 1, p = 1, 
J = 0.1 and q = 1). 



In the thin disk region, we may approximate cr r, while in the narrow jet region, 
(T ~ z. The system is mirror symmetry with respect to the z = plane, and the axes 
r = and z = are left as singularities. We consider -0 such that 



ip = iP{t) 



(21) 



where J and p {D and q) are positive constants, which control the strength of the jet 
(disk) flow. 

As shown in Fig. [51 this -0 models a disk-jet flow. The level sets of -0 (hence, those 
of r) are the streamlines of P . On the other hand, a serves as the coordinate directed 
parallel to the streamlines. We assume that pi is written as 

Pl(T,a) =pi(T)p||(a), (22) 

and, then, logpi = logp_L(T) + log p|| (ct). 

Let us see how the stream function defined by (j2ip satisfies the determining 
equations (fT5|). (fT7| . and ([T8| . i.e., we determine all other fields Pi, P2, v and 

h that allow this to be the solution. For arbitrary /(r) and (/(r), we observe 



V/-V.g = ^(T2 + l)/V. 



Hence, the left-had side of (|T7)) is (denoting g(r 

1 



L'^J — Vip ■ V log pi 



logp^(r)) 
[(r2 + l)V;" + 3ri/;'-(r2 + l)gV'] 



(23) 



For this quantity to balance with the right-hand side of pT)) . r~^pi must be a function 
of T if 12(0) ^ 0. Instead of demanding this relation for pi (cf. Rcmark[2|), we recourse 
to an assumption l2{ip) = (the implication of this simple condition will be discussed 
later). Then, (IT71) reduces to 

(t^ + 1)1/;" + StV'' - (r^ + l).g'7/'' = 0. (24) 
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Figure 3. The distribution of p± of the similarity solution (with D = 1, p = 1, 
J = 0.1 and 5=1). Contour levels are given in the log scale of p±. 



Wc note that the Beltrami condition is freed from p\\{(j)- This fact merits in 
solving ([151); sec ((27)) . 

For the specific form (|2ip of -0, we have to determine an appropriate g — logpj_ 
to satisfy i.e., 



9 =17 



ip" 3t 



i'' t2 + 1 

Jpip-l)TP+i + Dq{q + l) 3t 



Solving ([25]). we obtain 



and, thus. 



(25) 



9 logp± = log + - log(r + 1), 

{T^ + l)y^\JpTP+'^-Dq\ 
P± = ^5+1 ■ (26) 

In Fig. [21 we show the profile of p±{t). 

Remark 2 Here we considered the case of I'2{^j) = 0, hut a more general solution 
can he ohtained hy demanding r^'^pi to he a function of t. For example, we may put 
pi = P±iT)p\\ (cr) with p|| (tr) = = r^ + z"^ . Then, r~'^pi = (1 + t'^)pj_{t). 

5.2. Bernoulli relation in the disk region 

As mentioned above, this solution assumes l2{ip) (= A) = 0, and hence, I2 = P2rVg 
must uniformly distribute. In the disk region (the vicinity of z — 0), we may 
approximate Vg ~ \J MGjr (Keplerian velocity). Hence, p2 oc r^^/^. In Fig.[4l we 
show the profile of p = P\P2 for the case of p\ oc (i.e., p\\ ^constant). 
For p = P1P2 = p||/9^r^^/2, (ITSl) reads as 

.,_-d0/dz^_3/2^r-^ (27) 



2P||PJ- P\\{r) 
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along each streamline in the disk region. For a given j/, we can solve (|27p for pn to 
determine the density profile. 

In the disk region the Bernoulli relation ([T51) accounts as follows: by V • P = 0, 
we have Pr = pVr oc r~^. If P||(ct) = constant, for example, p oc 7*"^/^ (evaluated 
along a streamline in the disk region). Then, we have Vr — Vror~^^'^ with a (negative) 
constant Vro- Combining the azimuthal velocity Vg = Vgor^^^^ (which must be slightly 
smaller than the Keplerian velocity ^jMG/r) , we obtain 

By P2 oc r^^/^, we obtain dr{\ogp2) = ^(l/2)''^^- Hence, the Bernoulli relation p3| 
demands 

drh = {V^ ^ MG)r-^ , 

which yields h = {MG — V^)r~^. In this estimate, all components of the energy 
density (gravitational potential 0, kinetic energy V'^ /2 and enthalpy h) have a similar 
profile (oc r~^). 



5. 3. Bernoulli relation in the jet region 

In the jet region (vicinity of r = 0), the streamlines (contours of r = z/r) are almost 
vertical, and we may approximate a ~ z. 

Let us first estimate p2 using (|15p . which is approximated, in the jet region, by 

rpu = {logp2,'0} ~ (dri^) {dziogp2) 

= JpTP+^-d,{logp2), (28) 
z 

which shows that p2 is an increasing function of \z\. Using p — p|| (cr)pj^(T)p2, we 
integrate ([25)) along the streamline (t ^constant, cr « z): 

^ = -d f 1) = f r^+^ p^(r)p|| {z)z'dz. (29) 

With this ^2(2), we may estimate the toroidal (azimuthal) component of the velocity: 
Ve = l2/{rp2) ~ {l2T)/{zp2), where I2 and r are constant (the latter is constant along 
each streamline). We find that the kinetic energy V0/2 of the azimuthal velocity 
decreases as a function of \z\ (both by the geometric expansion factor z^^ and the 
viscosity effect p^^). The steep gradient of the corresponding hydrodynamic pressure 
yields a strong boost near the foot point (z « 0). 

The poloidal component of the kinetic energy is estimated as follows: We may 
approximate 



(30) 



2p2z4 ■ 

Here, the vertical distribution of the density p = p±{t)p\\ {<j)p2 is primarily dominated 
by p\\{(j) w p\\{z). 

At long distance from the origin, the jet has a natural similarity property. For 
simplicity, let us ignore the effect of viscosity (v = 0), and assume P2 = 1- Then, 
P|| oc \z\~^/'^ yields {V^ + V^)/2 oc z~^, which may balance with the gravitational 
potential energy (f) — —MG\z\~'^ . Note, that the azimuthal component of the kinetic 
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Figure 4. The contours (in log scale) of the density p in the similarity solution 
l|21|l of the Beltrami model (with parameters D = 1, p = 1, J = 0.1, and q = 1). 



energy disappears at large scale {Vg oc z~'^). The Bernoulli condition (fT3|) gives h 
that also has a similar distribution of oc 

Figure^] shows the profile of p = P±P\\P2 with pn oc (jet region) and 

P2 oc r^^/^ (disk region). Dividing P by the density p, we obtain the velocity field V; 
Fig. [1] shows the streamlines of V corresponding to Figs.[2]and|4l 

6. Summary and Concluding Remarks 

We have shown that the combination of a thin disk and narrowly-collimated jet is the 
unique structure that is amenable to the singularity of the Keplerian vorticity; the 
Beltrami condition — the alignment of flow and generalized vorticity — characterizes 
the geometry. Here the conventional vorticity is generalized as ^ to subtract the 
viscosity force causing the accretion and the centrifugal force of the Keplerian velocity. 
Identifying the disk-jet structure as a generalized Beltrami vortex, we will be able to 
understand the self-organization process in terms of the "generalized helicity." As we 
have learned in the present practice, the helicity of the generalized vorticity is the key 
parameter that characterizes the self-organizing of a disk-jet system. 

We end this paper with a short comment on the singularity of the Keplerian 
velocity. The similarity solution has a singularity at the origin (where 4> ~ —MGr~^ 
diverges), which disconnects the disk and jet parts of our solution. To "connect" both 
subsystems, we need a singular perturbation that avoids the divergence of physical 
quantities dictating the small-scale hierarchy on which the disk and jet regions are 
connected smoothly; Shiraishi et al. [22] describes how the different topologies of 
magnetic fleld-lines in the disk and jet regions can be connected in a "boundary layer" 
determined by the the Hall effect in a weakly ionized plasma. 
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Appendix A. MHD model 

While this article describes a pure fluid-mechanical model of jet collimation, many 
authors invoke a magnetic field, thrusting the center of the disk, to "guide" (and 
twist, as often observed) the flow of charged gas (plasma). Here the fluid vorticity plays 
the same role of a magnetic field. Indeed, the vorticity of the canonical momentum 
combines the fluid vorticity and the magnetic field: V x {mV + qA) [m is the mass, 
q is the charge, and A is the vector potential of electromagnetic field). Using this 
"canonical vorticity," we may readily extend the present model to include the effect 
of magnetic field (in a Keplerian system, however, the singularity of the fluid vorticity 
may be the principal part of the canonical vorticity) . 

The coupling of flow and magnetic field is described by the magnetohydrodynamic 
(MHD) equations: 

dtP + V ■ {VP) = {\7 X B)xB- p\/(j) - Vp - V ■ n, (A.l) 

dtB = V x{VxB), (A.2) 

where B is the magnetic field (the magnetic energy density \B\'^/8tt is normalized by 
the kinetic energy density Sq '■= PoK?/2). 

A stationary solution of (jA.2[) is given by 

B = ^P, (A.3) 

where /U is a certain scaler function (representing the reciprocal Alfvcn Much number). 
While (jA.Sp is not a general solution, other solutions are possible only if the electron 
pressure Pe or the electrostatic potential ip is huge (of the order of the kinetic energy 
density £q)] & perpendicular component of B with respect to P causes a Lorentz force 
on electrons, which must be balanced by a potential force V(pe — v)- If these energy 
densities are small, B of order unity is only possible in the parallel direction of P. 
Operating divergence on both sides of (|A.3[) . we find V • (/iP) = P ■ V/i = 0, implying 

Adding the magnetic field, the static force balance equation ([2]) is generalized as 

V ■ {VP) - [V X (/iP)] X (^P) = -pV0 - Vp - V • n. (A.4) 

The generalized vorticity is now combined with the magnetic field as 

1% := V X Pa - iip2V X B. (A.5) 

This generalized vorticity is compared with that of the Hall-MHD theory, i.e. the 
canonical vorticity |15| . Here the combination of the mechanical and electromagnetic 
components are scaled by the physical parameters p, (measuring the magnitude of the 
magnetic field) and p2 (reducing the mechanical component by the viscosity force). 
The Beltrami equations and (|17p are, respectively, generalized as 

Ap-ivV' = V(pr'/) -mV(Ai/), (A.6) 

(pr' - /^V2) + V^A • (Vpr' - mVp) = -Xpi^I. (A.7) 
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Combining (|A.6p and (IA.7I) . wc obtain 

7(1 -A^V) 

which reduces into (|17p when ^ = (i.e. unmagnetized) 
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